home *** CD-ROM | disk | FTP | other *** search
/ TPUG - Toronto PET Users Group / TPUG Users Group CD / TPUG Users Group CD.iso / COMAL / A-COMAL Series / (k)aac.d64 / basic-misshd3 (.txt) < prev    next >
Commodore BASIC  |  2007-02-28  |  3KB  |  152 lines
  1. 100 REM  2-D AIR-TO-GROUND MISSILE-FILE MISSHD3-BA  USE 1526/802 PRINTER
  2. 105 POKE53265,11: REM TURN 0FF SCREEN
  3. 110 OPEN3,4
  4. 120 CMD3
  5. 130 PRINT#3,CHR$(16)"082-D AGM-VERSION-5B 3/21/84-LINEAR-ATAU=.05,DEL BI=.0097"
  6. 140 PRINT#3,CHR$(10)
  7. 150 PRINT#3,CHR$(16)"02VA";
  8. 160 PRINT#3,CHR$(16)"17XT";
  9. 170 PRINT#3,CHR$(16)"31NL";
  10. 180 PRINT#3,CHR$(16)"46VT";
  11. 190 PRINT#3,CHR$(16)"61G2MT"
  12. 200 VA  = 1562
  13. 210 XT  = 5000
  14. 220 NL  = 4
  15. 230 VT  = 0
  16. 240 G2MT= 0
  17. 250 KA  =-2.125
  18. 260 PRINT VA,XT,NL,VT,G2MT
  19. 270 PRINT#3,VA,XT,NL,VT,G2MT
  20. 280 PRINT#3,CHR$(10)
  21. 290 PRINT#3,CHR$(16)"03T";
  22. 300 PRINT#3,CHR$(16)"09VM";
  23. 310 PRINT#3,CHR$(16)"17DX";
  24. 320 PRINT#3,CHR$(16)"25DZ";
  25. 330 PRINT#3,CHR$(16)"33GAMM";
  26. 340 PRINT#3,CHR$(16)"41THM";
  27. 350 PRINT#3,CHR$(16)"49NM";
  28. 360 PRINT#3,CHR$(16)"57ZD";
  29. 365 PRINT#3,CHR$(16)"65DEL"
  30. 369 PRINT#3
  31. 370 CLOSE3
  32. 380 DT=0.01
  33. 390 GC=32.1739
  34. 400 RAD=180/(null)
  35. 410 XM=0.
  36. 420 ZM=0.
  37. 430 ZT=500
  38. 440 SG=ATN((ZM-ZT)/(XT-XM))
  39. 450 SR=0
  40. 470 KR=5.
  41. 480 KT=5.
  42. 490 KN=KT*NL
  43. 500 G1MMD=0.
  44. 510 GAMM=0.
  45. 520 AL=0.
  46. 530 DEL=0
  47. 540 W=48.3
  48. 550 S=.196
  49. 560 CBAR=.5
  50. 570 IY=1.6
  51. 580 RHO=.0023769
  52. 590 CL=0.
  53. 600 VS=1116.9
  54. 610 VM=VA
  55. 620 G2MT=G2MT/RAD
  56. 630 G3MTD=0.
  57. 640 T5HDD=0
  58. 650 T6HD =0
  59. 660 THM  =0
  60. 670 T=0.
  61. 680 T2=.1
  62. 690 TR=385
  63. 692 ATAU=.05
  64. 694 BTAU=ATAU+DT
  65. 696 DEL=0
  66. 700 REM  START OF ITERATIVE CALCULATION
  67. 710 REM( TRACKING LOOP )
  68. 720 EP=SG-SR
  69. 730 S2GRD=KR*KT*EP
  70. 740 IF S2GRD >=.35 THEN S2GRD= .35
  71. 750 IF S2GRD <=-.35 THEN S2GRD=-.35
  72. 760 SR=SR+S2GRD*DT
  73. 770 D1EL=EP*KR*(KN/KA)-0.0097
  74. 775 DEL=(DEL*ATAU+D1EL*DT)/BTAU
  75. 780 IF DEL > .087 THEN DEL= .087
  76. 790 IF DEL <-.087 THEN DEL=-.087
  77. 810 REM( AERODYNAMICS AND CONTROL )
  78. 820 M=VM/VS
  79. 830 C1D=.669
  80. 840 C4DL=.045
  81. 850 C5MA=-3.85
  82. 860 C6MD=-3.52
  83. 870 C7MQ=-400
  84. 880 C2AL=13.6
  85. 890 C3LD=3.15
  86. 900 CL=AL*C2AL+DEL*C3LD
  87. 910 CD=C1D+C4DL*CL*CL
  88. 920 QA=.5*RHO*VM*VM
  89. 930 QS=QA*S
  90. 940 L=QS*CL
  91. 950 D=QS*CD
  92. 960 MTHD=QS*CBAR*CBAR*C7MQ*T6HD/(2*VM)
  93. 970 MAL =QS*CBAR*C5MA*AL
  94. 980 MDEL=QS*CBAR*C6MD*DEL
  95. 990 T5HDD=(MAL+MDEL+MTHD)/IY
  96. 1000 T6HD=T6HD+T5HDD*DT
  97. 1010 THM =THM+T6HD*DT
  98. 1020 G1MMD=((L/W)+(TR/W)*SIN(AL)-COS(GAMM))*GC/VM
  99. 1030 GAMM=GAMM+G1MMD*DT
  100. 1040 V1D=GC*((TR/W)*COS(AL)-D/W-SIN(GAMM))
  101. 1050 NM=VM*G1MMG/GC
  102. 1060 VM=VM+V1D*DT
  103. 1070 AL=THM-GAMM
  104. 1080 GS=SR-THM
  105. 1090 IF GS >= .35 THEN GS=0.35
  106. 1100 IF GS <=-.35 THEN GS=-.35
  107. 1110 REM( MISSILE-TARGET KINEMATICS )
  108. 1120 XM=XM+VM*COS(GAMM)*DT
  109. 1130 ZM=ZM-VM*SIN(GAMM)*DT
  110. 1140 G2MT=G2MT+G3MTD*DT
  111. 1150 XT=XT+VT*COS(G2MT)*DT
  112. 1160 ZT=ZT-VT*SIN(G2MT)*DT
  113. 1170 DX=XT-XM
  114. 1180 DZ=ZM-ZT
  115. 1190 SG=ATN(DZ/DX):IFDX<0THENSG=SG+(null)
  116. 1200 R=SQR(DX*DX+DZ*DZ)
  117. 1210 T=T+DT
  118. 1220 REM( LOGIC AND PRINTOUTS )
  119. 1230 IF R <= 30.GOTO 1250
  120. 1240 IF T <=(T2-DT) GOTO 1350
  121. 1242 A=VM*SIN(GAMM-SG)-VT*SIN(G2MT-SG)
  122. 1244 B=VM*COS(GAMM-SG)-VT*COS(G2MT-SG)
  123. 1246 DLT=R/B
  124. 1248 ZD=A*DLT
  125. 1250 PRINT INT(100*T)/100,INT(R)
  126. 1260 OPEN2,4,2
  127. 1262 OPEN1,4,1
  128. 1270 PRINT#2,"99.999"
  129. 1271 PRINT#1,T;
  130. 1272 PRINT#2,"99999.9-"
  131. 1273 PRINT#1,VM;
  132. 1274 PRINT#1,DX;
  133. 1275 PRINT#1,DZ;
  134. 1276 PRINT#2,"999.999-"
  135. 1278 PRINT#1,GAMM;
  136. 1280 PRINT#1,THM;
  137. 1283 PRINT#2,"999.99-"
  138. 1284 PRINT#1,NM;
  139. 1288 PRINT#1,ZD;
  140. 1290 PRINT#2,"99.9999-"
  141. 1292 PRINT#1,DEL
  142. 1300 CLOSE2:CLOSE1
  143. 1330 IF R <= 30 GOTO 1370
  144. 1340 T2=T2+.1
  145. 1350 IF T>= 5 THEN END
  146. 1360 GOTO 720
  147. 1370 OPEN3,4
  148. 1380 PRINT#3,"   ZD=";
  149. 1390 PRINT#3,INT(10000*ZD)/10000
  150. 1391 POKE53265,27: REM TURN SCREEN ON
  151. 1400 CLOSE3:END
  152.